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EVALUATION 


This  report  documents  an  electromagnetic  (EM)  fields  analysis  technique  that 
was  developed  for  a  class  of  bodies  that  can  be  modeled  as  rectangular  cylinders 
with  an  arbitrary  cross  sectional  geometry,  and  are  referred  to  as  "Bodies  of 
Translation"  (BOT).  This  technique  has  the  ability  to  model  EM  scattering, 
radiation  with  multiple  antennas,  aperture  coupling,  near  and  far  electric  and 
magnetic  fields  as  well  as  surface  current  distribution. 

This  technique  is  presently  limited  to  modeling  BOPs  which  have  open  ends 
and  slot  antennas.  A  subsequent  effort  will  address  these  limitations  and  develop 
procedures  for  modeling  the  ends  of  the  BOPs  and  off-surface  radiations  such  as 
monopoles  and  loops.  Emphasis  will  also  be  placed  on  techniques  which  will  allow 
one  to  hybrid  this  method  with  other  EM  analysis  techniques  in  order  to  model  the 
behavior  of  more  complex  structures  as  well  as  to  economize  computer  resources. 

DANIEL  E.  WARREN 
Project  Engineer 


1.  INTRODUCTION  AMD  BACKGROUND 


Modern  aircraft  and  missiles  require  complex  electromagnetic  (EM) 
systems  to  perform  their  roles  effectively.  To  characterize  the  behavior  of 
these  systems,  analytical  techniques  have  been  developed  to  predict  the  EM 
radiation  generated  and  scattered  by  aerospace  vehicles.  Recently,  predic¬ 
tive  techniques  based  on  the  method  of  moments  (MM)  have  been  developed  to 
solve  a  variety  of  antenna,  coupling,  and  field  penetration  problems.  In 
these  analyses,  the  radiating  or  scattering  structures  are  represented  by 
wires,  surface  patches,  and  wire  grids.  Because  of  the  computational  re¬ 
quirements  of  the  method  to  date,  application  of  these  methods  has  been 

2 

limited  to  bodies  with  surface  areas  on  the  order  of  X  .  Larger  surfaces 
2 

(of  'v*  45  X  )  can  be  analyzed  via  the  MM  technique  if  the  vehicle  body  has 

some  degree  of  symmetry ,  such  as  in  the  case  of  bodies  of  revolution  (BOR) . 

To  treat  complex-shaped  bodies,  such  as  parts  of  wing  sections  and  non- 

circular  aircraft  fuselages,  a  generalized  theoretical  formulation  called 

the  method  of  moments  for  bodies  of  translation,  abbreviated  here  as  MM/BOT, 
2  3 

has  been  developed.  *  This  formulation  treats  the  radiation  and  scattering 
from  bodies  with  active  and  passive  apertures.  In  addition,  the  technique 
has  been  extended  to  compute  the  fields  in  the  immediate  vicinity  of  a  slot 
antenna  as  well  as  fields  coupled  interiorly  through  rectangular  apertures 
in  the  BOT  surface.  The  MM/BOT  technique  combines  many  of  the  cost- 
effective  features  of  the  MM/BOR  analysis^  with  some  of  the  shape  flexibil¬ 
ity  of  the  wire-grid  approach  and  retains  the  ability  to  treat  difficult 
boundary  conditions  associated  with  realistic  radiating  and  scattering 
geometries. 


2.  SUMMARY  OF  COMPLETED  EFFORT 


The  principal  results  of  the  completed  effort  are  enumerated  below. 
Detailed  discussion  of  the  individual  topics  is  given  in  the  sections  indi¬ 
cated. 

o  The  MM/BOT  formulation  was  developed  for  a  generalized  BOT  configu¬ 
ration  to  treat  the  radiation  from  asymmetric  aperture  (slot) 
antennas  or  arrays  (Sections  4  and  5).  The  results  are  in  agreement 
with  the  predictions  of  the  MM/BOR  technique  for  asymmetric  slots 
embedded  in  a  BOR  surface  (Section  8). 

o  The  tM/BOT  formulation  was  extended  to  compute  all  the  electric  and 
magnetic  near-field  components  at  an  arbitrary  point  in  the  vicinity 
of  a  radiating  or  scattering  surface.  For  radiating  apertures,  an 
arbitrary  polarization  and  antenna  excitation  can  be  specified 
(Section  6). 

o  An  analysis  was  Implemented  to  determine  the  aperture-coupled  fields 
(both  electric  and  magnetic)  produced  by  EM  illumination  of  the  BOT 
from  an  arbitrary  angle  of  incidence  and  polarization  (Section  7). 

o  A  computer  algorithm  was  developed  to  Implement  all  aspects  of  the 
MM/BOT  formulation  in  a  hierarchy  of  user-oriented  computer  codes. 
The  codes,  written  in  FORTRAN  IV,  are  modular  and  machine  indepen¬ 
dent.  All  parts  of  the  codes  were  tested  and  installed  on  the  RADC 
computer  system.  A  user/systems  manual  (Volume  II  of  this  report) 
was  developed  for  the  MM/BOT  algorithm.  A  series  of  example  prob¬ 
lems  was  provided  to  illustrate  the  technique  for  the  prospective 
user. 

o  The  entire  computer  algorithm  was  tested  for  a  series  of  EM  problems 
and  validated  with  results  using  classical  boundary  value  solutions 
and  the  MM/BOR  codes  (Section  8) . 


3.  SUMMARY  OF  PREVIOPS  WORK 


The  original  formulation  of  the  method  of  moments  (MM)  was  applied 

initially  to  treat  the  radiation  and  scattering  from  thin  wire  structures*’^ 

4 

and  later  to  bodies  of  revolution.  As  shown  in  the  subsequent  discussion, 
the  MM  theory  can  be  extended  to  treat  finite  cylindrical  bodies  of  arbi¬ 
trary  cross  section,  denoted  here  as  bodies  of  translation  (BOT) .  Examples 
of  such  configurations  are  shown  in  Figure  1.  In  this  report,  the  cylinders 
are  assumed  to  be  uncapped  (i.e.,  open  at  the  ends).  The  theoretical 
development  for  this  class  of  bodies  parallels  in  part  the  MM/BOR  formula- 

4 

tion  and  retains  the  modal  expansion  concept  developed  in  the  latter  theory 

6  7  8 

Earlier,  Andreasen  ,  Wallenberg  and  Harrington  ,  and  Wilton  and  Mittra 

among  others,  have  treated  the  case  of  cylinders  of  arbitrary  cross  section 

but  of  Infinite  length.  The  case  of  finite-length  cylinders  of  arbitrary 

cross  sections  has  not  been  treated  previously.  Several  investigators 

examined  the  special  case  of  finite-length,  right-circular  cylinders  for 

9  10 

various  limiting  cases.  For  example,  Ufimtsev  ,  Kieburtz  ,  and 

Fialkovskii**  developed  solutions  for  thin  cylinders  with  ka  «  1;  while 
12 

Adey  considered  long  cylinders  with  ka  *  1,  where  ka  *  2TTa/X  and  a  is  the 

13 

cylinder  radius.  Williams  studied  the  diffraction  from  finite-length 

hollow  cylinders  where  the  open  ends  did  not  materially  influence  the  dif- 

14  15 

fracted  waves.  A  complete  study  of  tubular  cylinders  was  made  by  Kao  * 
for  arbitrary  length  and  kat  but  restricted  to  broadside  illumination.  His 
formulation  resulted  in  a  pair  of  decoupled  integral  equations  that  were 
solved  to  yield  the  axial  and  the  circumferential  currents  on  the  cylinder. 
Recently,  Davis  and  Mittra^  examined  the  current  distribution  on  an  open 
cylinder  of  1  X  length  and  ka  -  1  using  a  hybrid  formulation  incorporating 
the  electric  and  magnetic  field  integral  equation  representations  of 
Maxwellvs  equations  (i.e.,  EFIE  and  MFIE,  respectively).  The  present  MM/BOT 
formulation  treats  the  foregoing  problems  as  special  subcases. 


Figure  1.  Body  of  translation  (0OT)  configurations. 


4.  DEVELOPMENT  OF  MM/BOT  FORMULATION 

4.1  Electric  Field  Integral  Equation  for  BOT 

In  a  given  radiation  or  scattering  problem,  the  total  electric  field  Is 
given  as 


E  -  E1  +  E®,  (1) 

where  1  and  s  refer  to  the  Incident  and  scattered  fields,  respectively.  The 
scattered  field  in  turn  is  defined  in  terms  of  the  vector  potential.  A,  and 
scalar  potential $,  yielding 


E®  *  -  ju)A  -  V4>, 


(2) 


(3) 


(4) 


o 


(5) 


The  field  at  the  surface  of  the  conductor  can  be  expressed  in  terms  of  the 

^  S 

Incident  and  scattered  electric  fields  E  and  E  ,  respectively,  i.e., 


-  -  -  [jwA(J)  +  V*(J)]tan  =  L(J) , 


(6) 


where  J  is  the  unknown  current  density  on  the  surface  and  L(*)  is  a  linear 
integro-differential  operator  over  3.  Writing  out  Equation  (6)  explicitly 
yields 


(7) 
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4.2  Expansion  of  Surface  Currents 


Restricting  the  discussion  to  a  BOT,  Equation  (7)  Is  solved  for  J  by 
subdividing  the  domain  of  the  lntegro-dlfferentlal  equation  by  segmenting 
the  BOT  surface  S  into  length-wise  strips  as  shown  in  Figure  2,  On  the  con¬ 
ducting  surface  S  of  the  scatterer  (Figure  2)t  the  unknown  current  density  J 
is  expanded  in  a  double  sum  of  basis  functions  spanning  S,  described  by  the 
orthonormal  coordinates  (t,z),  so  that 

j .  p  ♦  ? .  y;  { v'jf .  <•> 

j  »n 


(*j.  Vj.  ijl  *  (tj.  Zj) 


(Xj.  Vj,  Zj)  =  (t,-,  Zj) 


Figure  2.  Segmented  BOT  configuration. 


Where  £  -  z/L,  x  -  t/T,  and  vn(5)  is  to  be  defined  later.  (L  ia  the  half- 
length  of  the  scattering  body  along  z,  and  T  is  the  girth  of  S  along  t 
divided  by  the  number  of  strips  chosen  to  represent  S.)  Equation  (8)  implies 
a  modal  expansion  of  J  along  £.  The  terms  f  ^  (x)  and  f * (x)  can  be  chosen  as 
pulse,  piecewise  sinusoidal,  or,  as  in  this  case,  triangle  functions,  l.e., 


1  -  Ix’Ux'l  <  1 

£?(X)  -  <  .  X»  -  T  -  X  ,  o  - 

J  (  0,  X*  >  1  j 


t  or  z. 


I 


For  a  nonclosed  BOX  (l.e. ,  the  silt  cylinder  In  Figure  1),  a  half-triangle 
function  Is  used  to  expand  the  currents  along  an  edge  parallel  to  z.  The 
choice  of  triangle  functions  allows  compact,  computer-coefficient  expres¬ 
sions  to  be  obtained  for  the  MM  Impedance  expressions.  Substituting  Equation 
(8)  into  Equation  (7)  and  forming  the  Inner  products  via  the  Galerkin  pro¬ 
cedure  with  respect  to  a  set  of  trial  functions  with  a  -  t  or  z, 

yields  the  matrix  equation. 


(10) 


In  Equation  (10) ,  the  column  vectors  V  and  I  refer  to  the  equivalent  voltages 
and  currents  on  S,  respectively,  and  where  the  i-th  element  of  V  in  terms  of 
the  inner  product  for  the  n-th  mode  is  defined  as: 

C  ■  &  <*L-  **> 

(il) 

-  fda  E*  •  \f* .  «  /*ds  E*  •  u  f?(T)  v  (C)» 

J  tan  ni  J  tan  a  l  n 

S  S 

where  Va. ,  with  a  ■  t  or  z,  denotes  the  a-directed  voltage  on  the  BOT. 
nl 


4.3  Derivation  of  the  Impedance  Matrix,  ZgQj 

The  8tructure  of  Z^^  is  given  in  Figure  3,  where  each  Zm  Is  a  full 
submatrix.  The  foregoing  matrix  equation  bears  a  canonic  similarity  to  the 
MM/BOR  result.  However,  in  general,  the  expansion  functions  in  Equation  (8) 
are  not  orthonormal  with  respect  to  the  integral  operator  L(*)  over  the  sur¬ 
face  of  the  BOT  as  in  the  BOR  case.  Thus,  in  principle,  all  the  Z^  matrices 
are  present  in  Z„_,  since  there  is  no  modal  decoupling.  The  elements  of  the 
partitioned  submatrices  of  can  be  computed  from  the  general  expresaion: 

fc?)-  -//<•//-  [i™.  *^r(r  •«>'>) 

S  S'  *- 


4xR  * 


where 


(a,  3  ■  t  or  z)  and  R  *  v(x 


-  x’)2  +  (y  -  y')2  +  (*  -  *')2  . 
8)1  and  the  fact  that  -  J*. 


Using  the  expansion  for  J  [Equation  (8)]  and  the  fact  that  Wni  -  7\ 

(*  denotes  a  conjugate)  and  noting  that  for  any  vector  A,  »  T  3A^/3t)  + 


L“^(3A  /3C),  Equation  (8)  becomes: 
z 


*  ^Kn6  f d8  I  ds'  (cos(v  “  v,)  fi(T)  fj(T,) 

'  '  s  s' 

-  4-  fJ(T')  »*<C)  Vn(«  », 


k]L  -  £  /"•/  d*'  £J(T)  <w>  ’ 

'  r  Of 


vn(C)  *, 


1 


jitns  f  d.  f  d.'  (v*«')vn(«  -  i  i;(f)i„(o) 
s  s’  K 

•  f*(T)  f*(x’)  * 


(15) 


Where 


and 


ds 


dxd£;  n 


* 


_1_ 

4tt 


e-jK  V<$2p2  +  (?  -  C')Z 

VfiV  +  (C  -  5')Z 


p  -  ^  ^(x  -  x')2  +  (y  -  y')Z  . 


(16) 


and  V  and  vf  are  the  angles  between  the  t-curve  and  the  x-axis  at  points  x,y 

zt 

and  x’,y’  on  the  BOT,  respectively.  The  expression  for  (Z^)^  is  identical 
to  that  given  in  Equation  (14),  with  m  and  n,  and  i  and  j  interchanged.  For 
evaluation  of  Equations  (13-15),  the  functional  form  of  v  (C)  ■  exp(jnirc)  was 
chosen.  The  triangle  function  fj(*)»  [Equation  (9)],  and  its  derivative  f.j(*) 
are  approximated  by  four  pulses ,  denoted  as  and  T^,  respectively,  with 
p  *  1,  2,  3,  A  where 


T 

P 


U-  b  b  r}  - 


T 

P 


The  expressions  for  the  (i,j)-th  elements  of  the  submatrices  of  are 

obtained  by  carrying  out  the  surface  integrals  in  Equations  (13-15) •  The 
corresponding  triangle  functions  are  depicted  graphically  in  Figure  A.  Thus, 


-  J2Kn6  E  [iV  cos(V  G 


(17) 


i-th  triangle 


function 


Pulse  approximation 


(i  +  1)-th  triangle 
function 


j-tn  triangle 
function 


BOT  generating  curve 


Figure  4.  Pulse  approximation  for  triangle  functions  on  BOT  surface. 
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where  following  the  convention  in  Reference  4,  (1, j)  correspond  to  the 

Indexes  p,q,  respectively.  The  term  is  the  engle  between  the  t -curve  end 

the  x-axls  at  point  x  y_  on  the  BOT.  The  function  G _  is  the  Green's 

p  p  mn 

function  Integrated  over  ds  and  ds'(  where  these  integrations  can  be  carried 
out  quaal-analytlcally,  as  discussed  below. 


4.4  Evaluation  of  Gmn  Functions 

The  evaluation  of  Gjl[l  encompasses  two  surface  Integrals,  written 
explicitly  as: 
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where  the  x'  and  x  integrations  are  carried  out  over  the  p  and  q-th  subseg¬ 
ments  (strips)  associated  with  the  i-th  and  j-th  strips  (see  Figure  4). 
Letting  K  •  .5  “  C'»  the  two  Integrations  over  £  and  reduce  to 
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where 
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R  -  Vfi2  p2  +  £2  and  p2  -  ~  [(x  -  x1)2  +  (y  -  y’)2]. 

The  Integral  over  x’  using  the  centroid  approximation  over  the  interval  Ax^’ 
yields 
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The  integration  over  T  can  be  carried  out  analytically  using  a  Taylor  expan 
sion  of  the  integrand  about  Rpq  and  completing  the  square  in  the  integrand. 
Using  only  one  term  in  this  expansion  yields 
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The  approximation  in  Equation  (25)  Implies  the  condition  that 


The  Integration  over  s  in  Equation  (25)  Is  carried  out  analytically.  Then 
approximating  the  integral  over  £  in  Equation  (24)  by  a  series. 
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where 
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For  the  self-terms  (i.e.,  p  »  0),  Equation  (24)  reduces  to  the  form, 
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The  Integral  over  £  can  be  evaluated  numerically  via  Gaussian  quadrature  or 
a  Simpson  integration  routine  after  an  integration  by  parts  to  soften  the 
In  £  singularity  at  the  lower  limit.  Thus, 
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This  completes  the  evaluation  of  the  elements  of  the  ZD/VP  matrix.  The 
analytical  and  computational  complexity  involved  in  this  analysis  is  approxi¬ 
mately  equivalent  to  that  encountered  in  the  MM/BOR  formulation.  Thus,  the 
matrix  fill  times  for  each  impedance  element  in  MM/BOR  and  MM/BOT  are  about 

the  same.  Furthermore,  detailed  examination  of  the  G  (i,i)  function  shows 

mn 

that  it  is  maximum  when  m  ■  n  and  i  *  j,  i.e.,  the  largest  values  occur  on 
the  main  diagonal  of  the  Zfi0T  matrix,  and  the  self-terms  contribute  the  most. 
These  properties  lead  to  a  diagonally  strong  overall  matrix  that  is  inter- 
coupled  significantly  only  for  the  neighboring  modes.  Computer  simulations 
.have  confirmed  this  feature. 


4.5  General  Structure  of  the  Impedance  Matrix 

Using  the  results  obtained  for  the  individual  members  of  the  partitioned 

submatrlce8  in  Equations  (17-20) ,  the  overall  impedance  matrix  for  a  BOT  can 

be  constructed  as  shown  in  Figure  3.  In  general,  this  matrix  is  full  since 

there  is  no  modal  decoupling,  (In  the  MM/BOR  analysis,  only  the  matrices  Z^ 

along  the  principal  diagonal  are  present.)  However,  certain  symmetries  exist 

for  the  G  and  Z  that  reduce  the  fill  time  of  the  Z^-  matrix.  Specifi- 
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where  . (G _ )  denotes  the  transpose  of  G  .  In  the  implementation  of  the 
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analysis,  only  the  lower  triangular  quadrant  of  the  Z BQT  matrix  is  computed, 
lee,,  the  partitioned  submatrices  are  filled  for  0  £  m  <  NMODE,  and 
-m  <*  n  m,  where  NMODE  is  the  total  number  of  axial  modes  (including  m  ■  0) 
used  in  this  analysis.  To  compute  the  inverse  of  the  entire  ZRQT 

matrix  must  be  filled,  which  can  be  accomplished  from  the  following  symmetry 
relationships: 


These  relations  are  exploited  in  the  MM/BOT  computer  algorithm  to  minimize 
the  fill-time  and  aid  in  the  solution  of  the  matrix  equations  arising  in  the 
formulation.  If  the  BOT  has  several  (physical)  planes  of  symmetry,  addi¬ 
tional  relationships  can  be  established  within  each  of  the  partitioned  sub- 
matrices,  again  allowing  savings  to  be  made  in  the  computational  requirements. 
In  summary,  the  MM/BOT  formulation  yields  an  overall  network  representation 
composed  of  diagonally  strong  matrices,  possessing  certain  symmetries  and 
appearing  to  be  lntercoupled  significantly  only  for  neighboring  modes. 


5.  PAR-FIELD  RADIATION  AND  SCATTERING  ANALYSIS 


Having  obtained  the  expressions  for  the  ZR0T  matrix,  the  currents  on 
the  body  can  be  obtained  by  solving  the  network  matrix  equation  for  the 
current  vector  I .  In  turn,  knowing  I,  the  radiated  or  scattered  fields  can 
be  obtained  as  in  the  MM/BOR  analysis  if  the  radiation  transfer  matrices 
Rt®,  RZ®,  Rt<^>,  and  RZ<^  are  given.  Formally,  the  i-th  element  of  the  trans- 
fer  matrices  is  defined  as 


<Oi  -<£i>  *?>  <“> 

where  the  superscript  a  denotes  t  or  zt  and  u  denotes  the  0  or  <p  polariza¬ 
tion  of  the  radiated  field  E^.  In  spherical  coordinates, 
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where  the  field  point  of  measurement  is  at  6^,  <J>r,  and  p  is  the  distance  to 
a  point  on  the  BOT  surface,  measured  from  the  origin  (see  Figure  5).  Expres¬ 
sing  the  Inner  product  in  Equation  (32)  explicitly,  the  transfer  matrices  are 
given  as: 
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Referring  to  Figure  5,  the  unit  vectors  for  the  BOT  geometry  are  given  by: 
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Using  a  centroid  approximation  for  the  t-integratlon,  a  pulse  approximation 

for  the  triangle  functions  f  .(•)  and  letting  v  (C)  *  exp(jnirC)  in  Equations 

x  n 

(34-36)  f  the  expressions  for  the  i-th  element  of  the  transfer  matrices 
become 
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where 
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Using  the  R-matrices  above,  the  far-fleld  power  radiation  patterns  and  the 
scattering  cross  sections  can  be  computed* 

5*1  Far  Fields 

The  total  radiated  far-field  in  u-polarization  can  be  computed  from 

—1  kr 
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where  the  sum  Is  taken  over  all  modes  m,n  used  in  the  analysis*  In  general, 
the  far-fleld  power  radiated  by  a  EOT,  excited  by  an  arbitrary  antenna  con¬ 
figuration,  is  given  by 


E  tf>B,„nvn]  . 
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where 
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where  the  caret  Indicates  a  row  vector.  The  matrix  Y  denotes  the  parti* 
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tinned  submatrix  corresponding  to  the  (m,n)-th  modes  in  the  Inverted 
matrix.  The  form  of  the  excitation  voltage  vector  depends  upon  the  type 
and  location  of  the  antenna  on  the  BOT.  As  an  example,  consider  a  series  of 
K  aperture  (slot)  antennas  embedded  in  the  BOT  surface  and  centered  at 
(t^  k  ■  1,  2  .  .  .  K.  Asstiming  that  the  apertures  are  rectangular, 

then  if  the  aperture  at  the  i-th  axial  strip  subtends  one  triangle  function 
f^x)  and  an  axial  width  of  -  Cq)  , 

v“i  -  TC  <^i-  *i> 

m  J  dT  f±(x)  y*  dc  uQ  •  E^T.C)  e'^n1T?,  a  -  t  or  z  (A5) 


In  general,  the  aperture  excitation  function  E^(t,£)  can  be  specified  to  be 
of  any  form.  For  this  discussion,  let  the  slot  be  uniformly  excited  by 
E^(V/m)  in  the  a-polarization,  so  that 
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In  the  above  expressions,  it  Is  assumed  that  the  1-th  slot  Is  excited 
uniformly  in  the  t  and/or  z  polarization  by  an  electric  field,  E^,  which  may 
be  a  complex  quantity.  In  the  program  if  5*  0,  then  (Ua)  ^  is  represented 
by  an  array  of  the  form 


(50) 


If  the  slots  are  excited  solely  in  the  t -polarization,  then  ■  0.  For 
simplicity,  U*  and  can  be  set  equal  to  unity.  By  obvious  extension,  these 
results  can  be  generalized  for  a  BOT  containing  a  series  of  slots  at  a  given 
axial  strip,  as  in  the  case  of  two-dimensional  antenna  arrays. 


5.2  Scattered  Fields 

In  general,  the  cross  section  in  terms  of  the  scattering  geometry  shown 
in  Figure  6  Is  given  as 


where  q  and  p  denote  the  polarization  of  the  incident  and  scattered  fields; 

(9. ,  <J>.)  and  (0  ,  <J>  )  are  the  incident  and  scattering  angles,  and  the  trans- 
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fer  matrices  Ra^  (a  »  t,z;  0  »  p,q)  express  the  relationships  between  the 
m 

current  on  S  and  an  observation  point  in  free  space.  Examples  of  cross 
sections  computed  from  this  expression  are  given  in  Section  8. 


Figure  6.  Coordinate  geometry  for  scattering  analysis. 


6.  NEAR-FIELD  ANALYSIS 


In  the  foregoing  sections,  the  radiated  and  scattered  fields  were  com¬ 
puted  at  field  points  sufficiently  far  from  the  BOT  so  that  the  wave  fronts 
were  planar  and  the  field  components  transverse.  Thus,  the  magnetic  fields 
could  be  obtained  from  the  computed  electric  fields  via  the  free-space  wave 
impedance,  n  “  377  ft.  In  this  section,  the  MM/BOT  formulation  Is  extended 
to  permit  the  electric  and  magnetic  fields  to  be  computed  at  points  <  X 
distant  from  the  BOT  surface.  First,  the  near-fleld  expressions  for  the 
electric  field  will  be  derived,  followed  by  the  corresponding  results  for 
the  magnetic  field. 

6.1  E-Near-Field  Formulation 

In  general,  the  electric  field  at  a  free-space  point  r',  resulting  from 
a  surface  current  density  J  on  a  surface  S,  is  given  by 

E(r')  ■  L(J) ,  (52) 

where  L(*)  Is  the  integro-dlfferential  operator  defined  in  Equation  (7).  In 
this  discussion,  the  current  density  can  be  induced  by  an  Incident  wave 
illuminating  the  body  as  in  a  scattering  problem  or  by  an  active  aperture  on 
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S  as  in  a  radiation  problem.  The  second  term  of  L(J)  in  Equation  (7)  can  be 
rewritten  using  the  relationship  that 
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where  4^(r  -  rf)  ■  [(1/R^)  +  j(k/R)]  #(r  -  r')  and  ♦  (r  -r')  is  the  free-space 

Green's  function.  If  S  corresponds  to  a  BOT,  the  electric  field  at  r'  can 
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be  rewritten  in  terms  of  the  t  and  z  directed  components  of  J.  Using  the 
modal  expansion  for  3  [Equation  (8)]  In  Equation  (52), 
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B(r  -  r')  ■  u  (x  -  x')  +  u (y  -  y')  +  u  (z  -  s'),  end  R  is  the  distance 
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between  field  and  source  points  defined  earlier.  The  evaluation  of  the 
integrals  in  the  foregoing  expression  follows  the  schema  used  in  determining 
the  integrals  in  the  impedance  expressions  (see  Equations  (13-15)].  Using  a 
pulse  approximation  for  the  f ^ ( • )  triangle  functions,  two  generic  integrals 
result  with  integrands  containing  4>(r  -  r*)  and  ^(r  -  r').  These  are  eval¬ 
uated  next.  Denoting  the  field  point  r'  with  (x^,  ?'), 


I,  -  f  dc  /*dx  e^nTr!*  $(r  -  r') 


where  R  »  Vfi2  p2  +  £2  with  6  and  p  defined  as  before  (Section  4.4).  The 
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T  Integration  is  carried  out  analytically  after  the  Integrand  Is  expanded  in 
a  Taylor  series  about  R^.  Then,  1^  becomes 
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and  f(£  )  and  its  associated  parameters  are  as  defined  before  in  Equations 
(27-28).  Note  the  evaluation  of  1^  is  similar  to  that  of  Gmw(p,q) ,  except  in 
the  present  case  there  is  only  one  surface  integral  to  evaluate  instead  of 


The  evaluation  of  the  Integrals  containing  ^(r  -  r')  parallels  the 
steps  followed  above.  Specifically, 


rift t  m  *+*»  -*>«>- 


I2(Y)  -  /  dc  / 
-1  j 


dt  e3™1’  (c  -  c')  *x(r  -  r') 


.  JnirC 


4  1-C  -jKR 


,  Y  ■  0  or  1, 


where  5  ■  C  -  Expanding  the  integrand  again  about  R^,  the  T  integration 
can  be  performed  as  before.  The  resulting  expression  for  I2(y)  becomes 
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where  X  and  all  the  other  parameters  have  been  defined  before.  The  expres¬ 
sions  for  and  have  no  self -terms  since  the  point  of  observation 
*f  (■  r  )  Is  assumed  to  be  different  than  r,  l.e. .  R  +  0. 


Using  the  results  for  1^  and  Ij ,  the  expression  for  the  electric  field 
In  tens  of  the  G  and  H  functions  become 
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6.2  H-Near-Field  Formulation 

In  general,  the  magnetic  field  is  given  in  terms  of  the  vector  poten¬ 
tial  A  and  the  incident  field  ^(rv)  as 

H(r’)  -  H1 (r')  +  V*  x  A(r’),  (68) 

Where  the  primed  coordinate  Is  at  an  arbitrary  test  point  where  the  field  Is 
to  be  sampled.  Letting  h  (r')  ■  0,  since  there  la  no  field  at  r't  except 
that  caused  by  currents  Induced  on  the  scattering  or  radiating  surface  S, 
Equation  (68)  becomes 


The  spherical  field  components  at  r'(-  p',  <p ' ,  6')  for  the  E  and  H  fields 
can  be  formed  by  combining  the  cartesian  components  in  Equations  (67)  or 
(71).  Thus,  letting  A  denote  the  E  or  H  fields,  the  appropriate  spherical 
field  components  are  given  by: 

A0f  ■  cos0'  cos*'  +  Ay  co 80 *  sin$'  -  Ag  sin0'  (72 

V  "  *  Ax  mSa*'  +  Ay  c08*'  (73 

Api  ”  Ax  co8<j)'  +  Ay  sin<J>'.  (74 

In  the  foregoing  discus sion,  the  electric  and  magnetic  near-fields  were 
sampled  at  a  point.  An  alternate  approach  was  also  considered  in  which  the 
fields  were  sampled  and  averaged  over  a  rectangular  patch.  This  latter  for¬ 
mulation  was  a  generalization  of  that  given  in  Reference  (17) .  The  predic¬ 
tions  from  the  patch  and  point  near-field  analyses  were  compared  for  a 
number  of  BOT  configurations.  Numerical  simulation  showed  that  in  regions 
where  the  electromagnetic  wave  departs  significantly  from  being  planar  (i.e. 
in  the  vicinity  of  active  or  passive  apertures) f  the  point-sampled  fields 
were  more  accurate,  particularly  for  aperture-coupled  fields.  Actual  com¬ 
parison  of  these  formulations  is  detailed  in  Section  8. 


7.  APERTURE  COUPLED  FIELDS 


The  problem  of  electromagnetic  fields  penetrating  through  small  a per- 
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tures  has  been  treated  previously  using  the  Bethe  small  hole  theory. 

Recently,  the  MM  technique  has  been  applied  to  symmetric  and  asymmetric 
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apertures  In  BOR.  9  A  computer  program  has  been  developed  by  Schuman 

to  treat  the  former  class  of  problems.  Here,  the  MM  technique  Is  applied 
to  the  case  of  rectangular  asymmetric  apertures  embedded  In  the  BOT.  The 
analysis  uses  the  Schelkunoff  equivalence  theorem,  which  replaces  the  exter¬ 
nal  sources  Illuminating  a  body  with  apertures  with  an  equivalent  problem 
having  only  aperture  current  sources.  From  the  aperture  currents,  the  near 
fields  Inside  the  body  can  be  computed.  In  the  present  discussion,  the 
aperture  Is  assumed  to  lie  anywhere  on  the  BOT  surface;  however,  an  aperture 
near  the  ends  of  the  BOT  can  lead  to  anomalous  unphyslcal  results.  The 
aperture  edges  are  taken  to  lie  parallel  to  the  z  and  t  coordinates  of  the 
BOT,  a  restriction  that  can  be  relaxed  at  a  cost  of  greater  analytical  com¬ 
plexity.  The  subsequent  discussion  assumes  a  single  aperture.  The  extension 

of  the  analysis  to  several  apertures  is  straightforward.  The  Schelkunoff 
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equivalence  theorem  is  discussed  next  for  a  general  scattering  surface, 
followed  by  its  application  to  a  BOT. 

7.1  Schelkunoff  Equivalence  Theorem 

A  graphical  statement  of  this  theorem  is  shown  in  Figure  7.  The  orig¬ 
inal  problem  of  external  fields  £  and  H  illuminating  a  body  with  an  aperture 
is  depicted  in  Figure  7a.  The  internal  (aperture  coupled)  fields  are  ^ 

H2.  With  the  aperture  covered  (Figure  7b),  a  current  J  is  induced  in  the 

region  of  the  covered  aperture  as  a  result  of  the  external  fields  2?  and  3? 

*  o 

In  Figure  7c,  the  equivalent  current  in  the  aperture  region  is  shown  as  (-J ) 
because  of  the  composite  external  electric  and  magnetic  fields,  i.e., 

-  Ej)  and  -  H^),  respectively.  A  simple  superposition  of  the  problem 
depicted  in  Figures  7b  and  7c  yields  the  original  problem  in  Figure  7a.  The 
results  of  the  Schelkunoff  theorem  will  be  applied  to  a  BOT  geometry  in  the 
following  section. 


(A)  Original  problem 


(B)  Aperture  covered  (C)  Equivalent  aperture 

currents 
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Figure  7.  Schdkunoff's  equivalence  theorem. 


7.2  Expansion  of  the  Aperture  Voltage 

Assume  that  a  rectangular  aperture  is  centered  at  and  Ca  on  the  BOT 
(Figure  8).  (The  normalization  of  the  coordinates  adopted  in  Section  4.1  is 
again  used*)  Let  the  BOT  surface  be  illuminated  by  a  plane  wave  (^,  5^) 
eminating  from  an  angle  (0^,  Then  the  field  in  V/m  at  the  q-th  strip 

spanned  by  the  j-th  triangle  function  fj(t)  subtending  the  aperture  is 
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Aq  Xq<T>)  1^)  •  (75) 


where 
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and  A  and  A  are  the  t  and  z-polarized  field  components  induced  across  the 

q  q 

aperture  are  due  to  the  incident  fields.  These  components  will  be  determined 
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Figure  8.  BOT  geometry  for  aperture  analysis. 


in  terms  of  the  Illuminating  plane-wave  incident  on  the  entire  BOT  surface. 

In  Equation  (75),  U(-)  is  the  unit  pulse  function  and  X^(t),  (a  *  t  or  z)  is 
the  sampling  function  of  the  fields  in  the  aperture.  In  this  formulation, 
these  components  are  defined  to  be  pulse  functions  (Figure  9)  and  are  similar 
to  those  used  by  Schuman  in  the  MM/BOR  analysis  .  Note  the  x“(-)  functions 
are  modified  near  the  ends  of  the  aperture  to  approximate  the  edge  behavior 
of  the  t-  and  z-directed  currents.  If  the  aperture  subtends  N  strips,  then 
Equation  (75)  is  generalized  so  that  the  total  anerture  field  is 


Triangle  functions  1  1  I  I  I 


N 

i  -  £  Eq  .  (76) 


The  induced  voltage  at  the  aperture  [Equation  (75)]  can  be  expanded  in  terms 
of  the  expansion  functions  used  previously  on  the  BOT,  i.e.. 
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Equating  Equations  (75)  and  (77) 


The  limits  of  integration  x+  and  t  span  the  opening  of  the  aperture  in  the 
t  direction.  If  the  pulse  function  X®(*)  also  overlaps  the  (j  +  l)-th 
triangle  functions,  then  Equation  (80)  also  yields  F® 

For  a  given  illuminating  field  E,  in  either  0  or  $  polarization,  the 
voltage  vector  Vr  is  specified  so  that  the  currents  on  the  BOT  can  be  com¬ 
puted  on  the  surface  with  the  aperture  covered  over  (i.e..  Figure  7b)  from 
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7.3  Derivation  of  the  Aperture  Admittance 

Equation  (81)  yields  a  relationship  between  the  modal  current  components 
on  the  BOT  corresponding  to  the  i-th  triangle  function  (i.e.,  I  n  ■  0, 

+1,  ...)  because  of  a  voltage  sampled  by  the  j-th  triangle  function  (i.e., 
Vnj).  Writing  out  Equation  (81)  explicitly  with  only  the  i-th  rows  and 
j-th  colums  retained. 


•  t 


For  the  subsequent  discussion,  It  Is  convenient  to  adopt  the  following  short¬ 
hand  notation  for  Equation  (84),  The  m-th  row  of  Equation  (84)  can  be  written 
as: 


(85) 


where  denotes  the  individual  partitioned  submatrices  in  Equation  (84) 

The  admittance  matrix  in  Equation  (84)  is  a  full  matrix  with  all  m,n  modes 
included.  If  the  sampling  points  p  and  q  lie  within  the  aperture  boundary, 
then  the  above  admittance  matrix  of  the  apertureless  BOT  can  be  related  to 
the  admittance  of  the  aperture.  This  case  is  considered  next. 

Formally,  the  currents  Cap  in  the  region  corresponding  to  the  p-th 
sampling  function  Xp  In  the  covered  aperture  region  can  be  expressed  as 
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where  the  currents  C  are  expanded  In  terms  of  modal  components  and  triangle 

functions,  where  the  aperture  admittance  (Ya)  is  determined  in  terms  of 
/  \  t  z  '  ^ 

\  ii/mn’  an<*  w*iere  Aq  an<*  Aq  are  t*ie  t_  aQd  z-polarized  field  components 
across  the  aperture  (see  Equation  75) .  The  negative  sign  in  Equation  (86a) 
is  due  to  the  ScHelkunoff  theorem,  i.e.,  the  equivalent  aperture  currents  are 
the  negative  of  the  currents  induced  on  the  apertureless  body  in  the  vicinity 
of  the  aperture  resulting  from  the  illuminating  fields  (i.e..  Figure  7c). 
Substituting  Equation  (85)  into  Equation  (86a)  and  using  the  results  of 
Equation  (83) , 
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where 


«.* 


is  the  (k,l)-th  element  of  the  entire  inverted  Z0/Vp  matrix  for 

BOX 


mode  pair  (m,n),  and  a, 8  denotes  tt,  tz,  zt,  or  zz.  The  bracketed  terms  in 

Equation  (87a)  constitute  a  matrix  partitioned  like  Y  in  Equation  (86).  The 

& 

expression  for  the  Individual  elements  of  the  admittance  matrix  is  obtained 
by  equating  the  right  sides  of  Equations  (86)  and  (87a),  i.e.. 
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Where  Dn  “d  F  are  *iven  in  Equations  (79)  and  (80),  respectively.  Thus, 
the  individual  elements  of  the  aperture  admittance  are  composed  of  a  triplet 
of  matrix  elements  from  the  inverted  Zfi0T  modified  by  the  integral  of  the 
aperture  sampling  functions,  represented  by  the  F  functions.  The  above 
expression  is  the  admittance  of  an  asymmetric  aperture  and  is  generically 
similar  to  the  BOR  results  given  in  Reference  22  for  the  circumferentially 
symmetric  aperture  problem. 


7.4  Equivalent  Aperture  Excitation  Voltage 


To  compute  the  aperture  coupled  fields  (i.e.,  E2  and  H2  in  Figure  7c), 
the  voltage  induced  in  the  aperture  as  a  result  of  the  illuminating  fields 
must  be  computed.  This  voltage  is  the  equivalent  aperture  excitation  voltage 
which  can  be  obtained  from  the  aperture  admittance  and  the  currents  in  the 
aperture  region.  For  the  p-th  current  C®  (a  -  t  or  z),  the  equivalent 
aperture  voltage  EV*,  sampled  by  X®(* )  is  given  by 
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where  K  denotes  the  number  of  pulse-sampling  functions  Xp(*)  used  to  describe 

the  current  in  the  aperture  region.  Having  determined  EV*  (q  ■  1,  .  .  .  K), 

01  ** 

the  voltage  excitation  V  .  corresponding  to  the  j-th  triangle  function  on  the 
original  BOT  spanning  Xp(- )  can  be  obtained  from 


'ev*  (F*  .  +  F*  +  Ffc  ,  )' 
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Then  the  effective  currents  on  the  BOT  in  the  presence  of  the  aperture  for 
the  m-th  mode  and  the  i-th  triangle  function  are  given  by: 


",J  .(Z«L  !  (z“). 


where  the  indexes  m  and  j  run  over  all  the  modes  and  triangle  functions  on 
the  BOT,  respectively.  The  resulting  total  currents  on  the  BOT  with  an 
aperture  is  reconstructed  from  the  modal  coefficients  so  that  at  the  i-th 
sample  point  the  current  is 


^(T)  eJ 


The  axial  distribution  of  the  currents  is  obtained  by  evaluating  Equation 
(92)  for  |e|  <  1.  At  the  edge  |c|  *  1,  theoretically,  an  infinite  number 
of  modes  is  required  to  obtain  an  accurate  representation  of  the  current. 
Since  only  a  finite  number  of  modes  can  be  used  in  practical  application 
of  the  MM/BOT  formulation,  care  must  be  taken  to  interpret  the  current  dis 
tributlon  near  an  edge. 


7.5  Computation  of  Aperture-Coupled  Fields 


in  Section  6.  The  aperture-coupled  electric  and  magnetic  fields  are  given 

by  Equations  (67)  and  (71),  respectively,  where  the  current  coefficients  1° 

u; 

are  replaced  by  given  by  Equation  (91) »  and  the  observation  point  r1  is 
within  the  BOT •  The  three  components  P,  0  for  the  electric  and  magnetic 
fields  are  again  obtained  from  Equations  (72-74). 


8.  VALIDATION  OF  THE  MM/BOT  FORMULATION 

The  MM/BOT  analysis  was  applied  to  compute  the  far  fields  radiated  and 
scattered  by  a  BOT  as  well  as  the  near  fields  and  aperture-coupled  fields 
for  various  antenna,  aperture,  and  body  configurations.  The  results  were 
compared  with  data  obtained  using  accepted  experimental  or  other  theoretical 
methods.  In  the  subsequent  examples,  the  results  for  the  scattering  cross 
sections,  near  fields,  and  aperture-coupled  fields  are  not  normalized. 


8.1  Validation  of  Far-Field  Analysis 

Because  of  the  availability  of  computer  codes  for  the  MM/BOR  formula¬ 
tion,^  the  far-field  radiation  patterns  were  computed  for  antennas  embedded 
in  a  BOR.  The  results  from  the  MM/BOR  and  MM/BOT  formulations  were  compared. 
As  an  example,  a  slot  antenna  embedded  in  a  right  cylinder  of  radius  0.216  X 
and  length  2.76  X  is  depicted  in  Figure  10.  The  ^-polarized  slot  subtends  an 
angle  of  45°  and  is  2.06  X  long.  The  pitch  (vertical)  and  roll  plane  power 
radiation  patterns  are  plotted  in  linear  power,  normalized  to  the  MM/BOR 
results.  There  is  excellent  agreement  between  the  MM/BOR  and  MM/BOT  results, 
with  seven  circumferential  and  four  axial  modes  being  used  in  the  respective 
calculations.  The  sensitivity  of  the  MM/BOT  calculated  pitch  and  roll  plane 
patterns  for  the  above  slotted  cylinder  as  a  function  of  modal  sparsing  is 
shown  in  Figures  11  and  12,  respectively.  The  calculations  show  that  use  of 
only  diagonal  modes  (l.e.,  m  *  n)  in  the  MM/BOT  results  in  an  approximate  10X 
deviation  from  the  MM/BOR  patterns.  The  power  distribution,  normalized  to  an 
isotropic  radiator,  versus  mode  number  is  shown  in  Figure  13  for  this  problem. 
(The  numbers  in  parentheses  are  exponents.  "Negative"  powers  arise  from  com¬ 
puter  round-off  errors  and  are  insignificantly  small.  Negative  powers  are 
often  also  obtained  for  certain  modes  in  the  MM/BOR  analysis  of  Mautz  and 
Harrington.)  The  maximum  powers  occur  in  the  self-modes  (i.e.,  m  ■  n) . 
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Figure  12.  Sensitivity  of  radiation  patterns  to  modal  sparsing: 


Slot  dimensions: 
2.06X  x  0.084X 
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roll  plane  (d-fod  axial  slot). 


The  applicability  of  the  present  analysis  is  to  degenerate  BOT  surfaces 
is  demonstrated  in  Figures  16  and  17,  where  the  monostatic  scattering  cross 
sections  for  a  flat  plate  (2  X  on  a  side)  and  the  bistatic  cross  section  for 
a  parabolic  cylinder  of  2,76  X  length  are  depicted.  In  both  cases,  the  sur¬ 
faces  are  illuminated  broadside  (i.e.,  <J>^  *  0,  0^  ■  90).  For  these  calcula¬ 
tions,  4  modes  and  16  triangle  functions  were  used.  The  MM/BOT  results  in 
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Figure  16  are  in  good  agreement  with  the  experimental  results  of  Ross. 
Similarly,  the  TM  polarized  results  for  the  finite-length  parabolic  cylinder 
in  Figure  17  are  in  good  agreement  with  the  analytical  predictions  of 
Andreasen. 


Broadside  0S  <deg)  End-on 

Figure  14.  Comparison  of  MM/BOR  and  MM/BOT  computed  bistatic  cross  section  for  open  cylinder. 
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Figure  15.  Bistatic  scattering  calculations  for  a  square  cylinder. 


Figure  16.  Monostatic  scattering  cross  section  for  a  square  plate. 


8»2  Validation  of  Near  Field  and  Aperture  Analysis 

Validation  of  the  near-field  formulation  for  the  electric  and  magnetic 
fields  is  shown  in  the  subsequent  examples.  In  this  phase  of  the  investiga¬ 
tion,  the  fields  were  computed  by  the  MM/BOT  technique  for  points  at 
distances  «  X  from  the  radiating  or  scattering  surfaces.  As  a  benchmark, 
the  exact  boundary  value  solution  for  an  infinite  right-circular  cylinder  fed 
with  a  ^-polarized  slit  was  used  (Figure  18). 27  The  slit  subtended  an  angle 
of  <pQ  degrees  and  was  excited  with  a  uniform  field.  If  the  field  is  sampled 
near  the  cylinder,  the  predictions  of  the  MM/BOT  analysis  for  a  finite 
cylinder  can  be  compared  with  the  exact  (classical)  solution  obtained  for 
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Figure  IS.  Classical  solution  for  slit  cylinder  *  excited  slit). 


the  Infinite  cylinder.  A  comparison  of  these  solutions  for  a  cylinder  with 
ka  -  1.35  and  a  slit  of  45°  is  shown  in  Figure  19.  For  the  MM/BOT  analysis, 
the  cylinder  length  was  2.76  X  long  and  the  slit  length  was  2  X.  The  near 
fields  are  sampled  along  a  line  bisecting  the  aperture,  resulting  in  ■  0. 
The  closest  field  sampling  was  at  1.2a  which  corresponds  to  0.2a  (-  0.04  X) 
from  the  plane  of  the  aperture.  The  calculations  using  four  and  seven  modes 
produced  practically  the  same  results.  Also  shown  in  Figure  19  are  the 
fields  obtained  from  a  patch  near-field  formulation  in  which  the  sampled 
fields  are  averaged  over  a  flat  strip.  The  large  discrepancy  of  the  patch 
results  from  the  exact  solution  at  points  near  the  body  is  due  to  the  fact 
that  in  this  region  the  EM  wave  departs  significantly  from  being  planar.  At 
distances  'U  10  X,  the  patch  and  point  formulations  coalesce.  (As  expected, 
at  these  distances,  the  classical  solution  for  the  infinite  cylinder  and  the 
B0T  results  for  the  finite  cylinder  diverge.) 

In  Figure  20,  the  three  near-field  components  sampled  at  a  radial  line 
at  <J>  -  45°  to  the  aperture  center  are  given.  Again,  the  results  of  the  exact 
and  the  MM/BOT  solutions  are  in  excellent  agreement.  The  corresponding 
results  for  field  points  sampled  at  $  *  90°  are  shown  in  Figure  21.  At 
sampling  distances  »  A  from  the  BOT,  decreases  and  the  EM  wave  front  tends 
to  approach  being  planar.  Finally,  the  near  fields  for  a  slit  subtending 
22.5°  are  shown  in  Figure  22.  In  these  calculations,  4  modes  were  used  and  33 
points  defined  the  circumference  of  the  cylinder.  Again  there  was  close 
agreement  between  the  exact  and  MM/BOT  solutions. 

An  application  of  the  aperture-coupled  field  analysis  of  Section  7  is 

depicted  in  Figure  23.  The  internal  fields,  sampled  along  a  radial  line 

bisecting  the  aperture,  are  induced  by  a  broadside  TM  Illumination  of  a  right 

circular  cylinder.  In  the  BOT  calculations,  the  cylinder  length  was  5.52  X, 
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with  the  aperture  subtending  22.5°  and  an  axial  length  of  4.96  X.  Senior 

considered  the  infinitely  long  cylinder  with  a  20°  infinite  slit  and  computed 

only  the  axial  electric  field,  E  •  As  can  be  seen,  the  BOT  analysis  is  in 

z 

good  agreement  with  this  result. 
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Figure  20.  Computation  of  near  fields  for  slit  cylinder  at  0  *  46°  (slit  angle  3  46°) 
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Figure  21.  Computation  of  near  fields  for  slit  cylinder  at  0  ■  90°  (slit  angle  *  45°) 
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Figure  22.  Computation  of  near  fields  for  slit  cylinder  at  0  -  90°  (slit  angle  •  22S°) 


The  MM/BOT  formulation  described  In  the  preceding  sections  was 
Implemented  with  a  computer  algorithm  which  is  described  in  detail  in 
Volume  II  of  this  report.  The  overall  structure  of  the  program  flow  compares 
with  that  of  the  MM/BOR  codes.  The  computational  complexity  is  approxi¬ 
mately  equivalent  for  both  MM/BOT  and  MM/BOR  for  a  given  size  body.  The 
matrix  fill-times  are  comparable.  The  major  differences  lie  in  the  fact 
that  the  modes  in  the  present  analysis  do  not  decouple  as  in  the  MM/BOR, 
although,  in  general,  the  resulting  network  matrices  remain  diagonally  strong 
and  have  certain  symmetries  [i.e..  Equations  (30-31)].  For  sufficient 
accuracy,  the  off-diagonal  submatrices  for  m  ^  n  can  sometimes  be  deleted 
from  the  computation  without  excessive  error  penalty  (i.e..  Figures  11-12). 

Adequate  computational  accuracy  is  achieved  when  the  BOT  surface  is 
segmented  into  strips  <  0.15  X  in  width.  The  number  of  axial  modes  chosen 
is  dependent  upon  the  spatial  accuracy  desired  for  the  surface  currents. 

While  the  examples  shown  in  the  validation  section  involved  mostly  right- 
circular  cylinders,  the  present  formalism  is  capable  of  treating  any  asym¬ 
metric  BOT,  such  as  a  wing  section.  An  example  of  this  case  is  given  in 
Volume  II. 
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